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Motivated by the general problem of extending the classical theory 
of holomorphic functions of a complex variable to the case of quater- 
nion functions, we give a notion of an W- derivative for functions of one 
quaternion variable. We show that the elementary quaternion func- 
tions introduced by Hamilton as well as the quaternion logarithm 
function possess such a derivative. We conclude by establishing rules 
0^ ' for calculating H-derivatives. 

^D , 2000 Mathematics Subject Classification. 30G35. 30A05, 30B10. 

l/j i Keywords. Quaternion, H-dcrivativc, elementary quaternion functions. 

en 

/vi ! 1 Introduction 

The importance of the theory of holomorphic (analytic) functions of one 
complex variable suggests looking for a similar theory for functions of three 
and more real variables. Considering the case of four variables, where the 
quaternion algebra H should replace the field C of complex numbers, is, of 
C3 . course, especially natural. 

Let us recall that the quaternion algebra was introduced by W.R. Hamil- 
ton in 1843 (see, e.g., 0H]), and that, according to Frobenius' theorem (see, 
e -g-) [13])) every finite-dimensional (associative) division algebra over the 
field K of real numbers is isomorphic either to R, or to C, or to H. 

In a sense, there are three well-known methods for constructing the 
theory of holomorphic functions of one complex variable: the derivative 
method, the polynomial method and the gradient method. In the case 
of quaternion functions, none of them leads to a satisfactory quaternionic 
analogue of the notion of a holomorphic function. 

Derivative method 

This method is based on the use of the limit definition of a derivative. In the 
case of quaternion functions, such an approach yields at least two notions 
of the quaternion derivative of a quaternion function f(z). The point is 
that, as quaternions do not commute, there are two different possibilities 
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for regarding the ratio ^ l~ : one C0lu d replace it either by [f(z + 
h) — f{z)]h^ 1 or by h^ 1 [f(z + h) — f(z)]. This leads to the notion of a 
right-hand derivative 

A(z) = lim[/(z + h)- /(z)}^ 1 

h— fO 

and to the notion of a left-hand derivative 

B(z) = }imh- l {f(z + h)-f(z)], 

provided that the corresponding one-sided limit exists. Actually, both of 
these notions are too restricted. It turns out [16j ESI HI [21j HH US] that 
only the functions ip(z) — az + b possess the right-hand derivative, and only 
the functions ip(z) = za + b possess the left-hand derivative, while only the 
functions x( z ) = rz + b possess both left- and right-hand derivatives, and 
they are equal. Here a and b are quaternions, while r is a real number. 

Polynomial method 

Let us consider a polynomial P(x, y)= Y^ A m . n x m y n of two real variables 

x, y with complex coefficients A m ^ n — u m . n + if3 m , n . Replacing x with 
i(z + z), and replacing y with ^ i(z — z), we obtain the polynomial P*(z, ~z) 
of the complex variables z = x + iy and z = x — iy. In order that P*(z,~z) 
be a polynomial of only the variable z, it is necessary and sufficient that 
P*(z,z) satisfies the well-known Cauchy-Riemann condition. This condi- 
tion is precisely what is needed for distinguishing the class of polynomials 
of the variable z. 

An analogous approach to the polynomials of the variables xo, x\, X2, 
X3 with quaternion coefficients is too general to give a desired result 28J. 
Indeed, the HausdorfF formulas xo — -Az — i\zi\ — i 2 zi 2 — izzi^), x\ = 
-^■(z - i\zi\ + i 2 zi 2 + hzi 3 ), x 2 = -^(z + i\Z%\ - i 2 zi 2 + izzh) and x 3 = 
jr-(z + i\zi\ + i 2 zi 2 — i^zis) 12! allow us to express the real coordinates Xk 
of the quaternion z — x$ + x\ii + x 2 i 2 + x^i^ in terms of z itself, without 
using the conjugate quaternion z — xq — x\i\ — x 2 i 2 — x-$iz- (Note that 
along with the HausdorfF formulas , we also have the formulas £0 = \ (z + z) , 
x\ — \[i\Z — zii), x 2 — \{i 2 ~z — zi 2 ) and x$ — \(i^z — zi^), but unlike 
the complex case, they are not essential.) Hence the functions which can 
be represented by quaternionic power series are just those which can be 
represented by power series in four real variables. 

Gradient method 

Loomann 18 and MenchofF 22 proved that any complex-valued continuous 
function that satisfies the Cauchy-Riemann condition in a complex domain, 
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is holomorphic in the same domainj [27j . Thus, in order to extend the the- 
ory of holomorphic functions to the quaternionic setting, one of the possi- 
ble ways is to try to find a quaternionic analogue of the Cauchy-Ricmann 
equations. In 1935, Fuetei|j[8 proposed such a quaternionic analogue by 
introducing two quaternion gradient operators 

d r d d d d 

■7T- = q h q H + q «2 + q »3 

OZ OXq OX\ OX2 OX3 

and 

d l d d d d 

tt = q hii-g M2-5 M3-5 — ■ 

az OTo asi aa;2 0x3 

He calls a quaternion function f(z) of the quaternion variable z — xq + 
xiii + X2^2 + £3*3 right-regular, provided 

d r f 

-37'°- < L1 > 

Left-regularity of / is defined in a similar way by requiring that 

d l f 

-87 = °- ^ 

f is called regular if it is simultaneously left and right regular. Such quater- 
nion functions are to be thought of as the appropriate generalization of 
holomorphic functions to the quaternionic setting. 

Any right or left regular quaternion function is harmonic, i.e. satisfies 
Laplace's equation 



d 2 d 2 d 2 a 2 

dx 2 , dx\ dx\ dx\ 



A^(z)=0, &= + + + 



Moreover, for any real harmonic function, there is a regular function whose 
real part is exactly the given function (see [Sj). Thus there are plenty of 
regular functions. 

Although Fueter's approach is quite powerful and gives rise to a fully 
formed theory of regular functions, it has some significant flaws. One such 
is the fact that even the functions ip n (z) = z n fail to be regular (they 
are non-harmonic, since, for example, A(z 2 ) = —4 and hence cannot be 
regular). Note that the functions Az n are regular [S]. Another flaw is that 
the class of regular functions of a quaternion variable do not form an algebra 



1 Tolstov 1 30! proved the some result replacing continuity of a function by its bound- 
edness. 

2 Dr. Rudolf Fueter is the author of the book "Synthetische Zahlentheorie" , Berlin 
and Leipzig, 1921, pp. VIII+271. 
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in the same sense that the holomorphic functions do: for example, regular 
functions cannot be multiplied to give further regular functions. 

The aim of the paper is to propose a new definition of a derivative 
for quaternion functions of one quaternion variable and show that all the 
elementary functions as well as the quaternion logarithm function possess 
such a derivative. We also give rules for calculating such derivatives. 

We conclude Introduction with an incomplete list of references where 
various applications of quaternions are discussed [21 [3j [6] [101 E2 H~5l 

El EH 123123 ESI ED- 

2 The notion of an H- derivative 

We begin by the following 

Definition 2.1. A quaternion function f{z), z = xq + x\i\ + X2I1 + X3J3, 
defined on some neighborhood Gci of a point z° = x { ^-\-x\ii+x t 2 ) i2 J rX S li^, 
is called H-differentiable at z° if there exist two sequences of quaternions 
Ak(z°) and Bk(z°) such that ^2 Af c (z a )Bk(z°) is finite and that the incre- 

k 

raent f(z + h) — f(z°) of the function f(z) can be represented as 

f(z° + h)- f(z°) = ]T A k (z a ) ■ h ■ B k {z a ) + u(z°, h), (2.1) 



k 



where 



^ ! ^f M =0 (2.2) 

h->0 \h\ 



and z° + h e G. In this case, the quaternion ^2Ak(z°)Bk(z a ) is called the 

k 
H -derivative of the function / at the point z° and is denoted f'(z). Thus 

f(z )=Y,Mz°)Bk(z ). (2.3) 

fc 

The uniqueness of the H-derivative follows from the fat that the right- 
hand part of (|2.3j) . if it exists, is just the partial derivative f' Xo {z°) of f(z) 
at z° with respect to its real variable. 

In the sequel, the symbol o{h) denotes any function u)(z , h) satisfying 
(I2~2l). 



Remark 2.2. Note that the same definition still makes perfectly good sense 
for any map between Banach algebras. Moreover, all the proofs of our results 
remain valid (except for Proposition 13. 4[ which still remain true if we take 
(p to be invertible in a neighborhood of z ), since they require only those 
properties of H which are also possessed by any Banach algebra. 
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The following functions introduced by Hamilton 
z n , n = 0,1,2,... , 









z 2 


z 3 


e* = 


= 1 + 


z + 
z 2 


2! + 

7 4 


3! 


sz = 


= 1 - 


2! 
z 3 


+ 4!- 
z 5 





8mz = z __ + __... 

are the 6asic elementary quaternion functions of one quaternionic variable z. 
Let us now show that the basic elementary functions are H-differentiable. 

Proposition 2.3. (z n )' = nz ll ~ l for n = 0, 1, 2, . . . , and for zeH. 

Proof. We first show that the following equality holds for n — 1,2,... 

(z + h) n -z n = z n - l h + z n - 2 hz + z n - 3 hz 2 

+ ■■■ + zhz n - 2 + hz n ~ l + o(h) (2.4) 

For n = 1 it is obvious. Assuming now that it is true for n = k, we find 

(z + h) k+1 - z k+1 = {z + h){z + h) k - z k+1 

= (z + h)(z k + z k - l h + z k - 2 hz + ■■■ + zhz k - 2 + hz k - x + o(h)) - z k+1 

= z k+1 + z k h + z k - 1 hz + ■■■ + z 2 hz k - 2 + zhz k - x + hz k + o(h)) - z k+1 

= z k h + z k - 1 hz + z k - 2 hz 2 + ■■■ + z 2 hz k - 2 + zhz*- 1 + hz k + o(h). 

It then follows from (f2~3]) and (|2~4]l that 

(z n )' = z 71 - 1 ■ 1 + z n - 2 ■ z + z n - 3 ■ z 2 + ■ ■ ■ + z • z n - 2 + 1 • z n - 1 = nz"- 1 . 

Thus by induction we have proved that (z n )' = nz 71 ^ 1 for all n = 
0,1,2,... Q.E.D. 

In order to proceed further, we need the following lemma. 
Lemma 2.4. The following equalities and estimates are valid for \h\ < 1: 

(z + h) 2 — z 2 zh + hz 



2! 2! 

(z + h) 3 - z 3 _ z 2 h + zhz + hz 2 

3! ~ 3! 

(z + h) 4 - z 4 _ z 3 h + z 2 hz + zhz 2 + hz 3 

4! = 4! 

(z + hf - z 5 _ z 4 h + z 3 hz + z 2 hz 2 + zhz 3 + hz 4 

5! = 5! 



.2 

A 4 



■A* 
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and so on, where 

A 2 = ^h 2 , \A 2 \<\h\ 2 , A 2 = o(h); 

A 3 = i (zh 2 + hzh + h 2 z + h 3 ), \A 3 \ < L {\z\ |/j| 2 + \h\ 3 ) 



2_\ h \2 {1+ ^ )< 2_ W 2._ 1 _^ f0J , N<1) 

f | Z ||/ l |2(l + |/,|)<|i| 2 ||/ l |2._^_ /or N > 1; 

1 



At = -|y (z 2 /i 2 + z/12/1 + zh, 2 z + zh 3 + hz 2 h + hzhz + hzh 2 + h 2 z 2 
+ h 2 zh + h 3 z + h 4 )), I Ail < |f (|z| 2 N 2 + |z| |^| 3 + H 4 ) 



< 



| r | /l | 2 (l + |/ l | + | /l |2 )< | r | /l |2. T ^_ [ /flr | z]<1) 

*\z\*\h\*(l + \h\ + \h\*)<£\z\*\h\*- T ^ for \z\ 



>1; 



A 5 = — (z 3 h 2 + z 2 hzh + z 2 h 2 z + z 2 h 3 + zhz 2 h + zhzhz + zhzh 2 

+ zh 2 z 2 + zh 2 zh + zh 3 z + zh 4 + hz 3 h + hz 2 hz + hz 2 h 2 
+ hzhz 2 + hzhzh + hzh 3 z + hzh 3 + h 2 z 3 + h 2 z 2 h + h 2 zhz 
+ h 2 zh 2 + h 3 z 2 + h 3 zh + h 4 z + h 5 ), 

\A 5 \<^(\z\ 3 \h\ 2 + \z\ 2 \h\ 3 + \z\\h\ 4 + \h\ 5 ) 

' ff \h\ 2 (l + \h\ + \h\ 2 + \h\ 3 ) < | \h\ 2 ■ j^ for \z\ < 1, 

2 ^\z\ 3 \h\ 2 (l + \h\ + \h\ 2 + \h\ 3 ) 

<^\z\ 3 \h\ 2 - J ± Kl for\z\>l, 



5! II II 1-|/»| 

and so on, 



(21 
A < I ™ ! 

IfNr 2 !^ 2 -™ for\z\>l- 



\h\ 2 -T^ for\z\<l, 



Therefore 



, N 2 • T^TTf • E h for \z\ < 1, 

I N 2 • t=W • E Ir M n " 2 V 1*1 > i 

n— 3 



and the series ^2 ^r ari ^ 2 ^r M™ 2 are converging by virtue of the ratio 

n— 3 n— 3 

oo 

£es£ j20j . Thus, ^2 \An\ = o(/i) for any fixed finite quaternion z. 

n— 3 
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Proposition 2.5. We have the equality 

(e z )' = e z . 

Proof. The equality 

z 2 z 3 
6 = 1 + Z+ 2! + 3! 
implies that, for any /i£i, 



z+h 



(z + h) 2 -z 2 (z + h) 3 -z 3 (z + h) A - z 4 



2! 3! 4! ' 

and applying Lemma 2.3 to the right-hand side of this equality, we obtain 

e z+h -e z = h + - ] {zh + hz) 

+ — (z 2 h + zhz + hz 2 ) 

+ — (z 3 h + z 2 hz + zhz 2 + hz 3 ) 
4! 



o(h). 



Therefore 

e z+h - e z = ( I + f { + ^ + ■ ■ ■ ) h 



hz+[^ + ^7 + ^--\ ) hz 2 



and hence 



\ z z 2 \ , / 1 

2! + 3! + 4! + '"r Z+ 3! ' 4! ' 5! 



(e z )' = l + 7TT + TT7 + 77 



o(h) (2.5) 



1 + 


z 
2! 


z 2 
+ 3! 


+ 


z 3 
4! 


+ 


2; 
2! 


z 2 
+ 3! 


+ 


z 3 
4! 






z 2 
+ 3! 


+ 
+ 


z 3 
4! 
z 3 
4! 


z 
2 2! 


z 2 
+ 3 3! +4 


z 3 
4! 


1 + 


2; 
1! 


z 2 
+ 2! 


+ 


z 3 
3! 



Q.E.D. 
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Proposition 2.6. The equality 

(sin z) = cos z 
is valid. 
Proof. 

sin(z + h) — s'mz 

. ,. (z + h) 3 {z + hf 
= + h) - v - ; + - — — z + 



= h- 



3! 5! 3! 5! 

(z + /i) 3 - z 3 (z + hf - z 5 
3! + 5! 



Therefore 



-h-^h + zhz + h**) 

+ — (z 4 /i + z 3 hz + z 2 hz 2 + zhz 3 + hz 4 ) -\ h o(h). 

5! 



z 4 



sin(z + /i)-sinz = /i+(-— + — )/i 



Hence 



*M-^ + ^]+M-ir + ^r)--+»w- 



z 4 



1 3! + 5! 3! + 5! 3! -1! 



1 "¥ + 4! 



Similarly, one has 
Proposition 2.7. TTie equality 

(cos z)' = — sinz 
is valid. 



Q.E.D. 
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3 Rules for calculating H-derivatives 

The rules for calculating H-derivatives are identical to those derived in a 
standard calculus course. 

Proposition 3.1. Let f and ip be two functions defined on a neighborhood 
of z° € H. If both f and ip are M-differentiable at z , then 

1. both cf and fc are M-differentiable at z° for all c € M and (cf)'(z ) = 
cf'(z°)and(fcy(z ) = f'(z )c; 

2. f + ip is M-differentiable at z° and (/ + tp)'(z°) = f'(z°) + tp'(z°); 

3. ftp is M-differentiable at z° and (ftp)'(z°) = f'(z°)ip(z ) + f (z)tp' (z°) . 

Proof. The proof of 1. is obvious. 

Since / and ip are H-differentiablc at z°, there are representations 

f(z° + h)- f(z°) = J2 AkhBk + °^)' 



y 



(z° +h)- p(z°) = J2 °khD k + o(h). 



Then 



(/ + <p)(z° +h)- (/ + <p)(z°) = [f(z° + h)- f(z )} + [p{z° + h)- ^(z°)] 

= Y, A k hB k + Y C k hD k + o{h), 



and hence 



(/ + p)'(«°) = J2 A kB k + Y, CkD k = f'(z°) + <f/(z°). 



This proves 2. 
Next, since 

f(z° + h)ip(z° + h)- f(z°Mz°) = 

= [f(z° + h)- f(z°)Mz° + h) + f(z°)[<p(z° + h)- <p(z )} 



Y A k hB k + o{h) 



p(z" + h) + f(z°) 



J2 C k hD k + o(h) 



J2 A khB k + o{h) ■ ^{z°)+J2 C khD k + o(h) 



f(z°) 



J2CkhD k + o(h) 



J2 A k hB k ) p(z°) + f(z°) J2 CkhD k + o(h), 
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it follows that 

(M(z°) = (5> fe B fc W°)+/(*°)£c fcJ D fe =f(z Mz )+f(z°)<p'(z ), 
^ k ' k 

proving 3. q.e.d. 

The following two corollaries are immediate. 

Corollary 3.2. 7f /i, /a, • • • ,/n ar ^ M-differentiable functions at a point 
Z , then their product /1/2 • • • fn * s also M-differentiable at z° and we have: 

(/1/2 • • • fn)'(z°) = f[(z°)h(z°) ■ • • f n (z°) + h{z°)f' 2 {z°)h{z°) ■ ■ ■ f n (z°)+ 

Corollary 3.3. If a function f is M-differentiable at a point z°, then f n is 
also M-differentiable at z° for all n = 1, 2, . . . and we have: 

cnV) - f'(z o )r- i (z°) + /(z°)/'(z°)r- 2 (z°) + . . . + r~ V)/'(A 

Proposition 3.4. If a function <p is M-differentiable at a point z° and if 
p 7^ in a neighborhood of z° , thei\j — is also M-differentiable at z° and 
we have: 

lV, n, 1 1 



ip ) <p(z°) (p{z") 



Proof. We first observe that for any two nonzero quaternions q\ and (72, the 
following equality 

Qi 1 - Q.2 1 = <h x {<h - <ta)<h x {Qi - ^q^ 1 - 9 2 " 1 (q , i - ffiOgf 1 

holds. Indeed, using that q^ is the inverse of q\ and q 2 is the inverse of 
q2, we obtain 

Qi l (qi - Q2)q2 1 (qi - ftOfe 1 - q^ l {q x - 92)^ 1 = 

= (1 - <h l( &) c h l { < li c h l - 1) - (fe V - l )l2 l 
= i<h l - <h l ){<li<h l ~ 1) - (feVife 1 - c h l ) 

= {l2 l 1\l2 l - <h X - 12 1 + Ql 1 ) - (feVfe 1 - C l2 1 )-. 

as needed. 



3 For each quaternion 9^0, there is a (unique) quaternion -, called the inverse of q, 
for which q ■ - = 1 = - ■ q. The inverse of q is sometimes denoted by the symbol q~ . 
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Putting tp(z + h) and <p(z°) in the equality, we obtain 
1 1 



<p(z° + h) <p(z°) 
1 



(p(z°) ip{z 
We now calculate 
1 1 



o^ W + h) - „(*»)] _^ y }. [ ^ +h )-^)] -1_ 



^o)^°^)-^°)]^o) 



+ R^W ^° + *> ~ ^ ^y [**° + *) - *(*°)] ^y 

= -^u)[E^H^) +0(/l) - 

Hence 

(£) (z ° ) = ~rW E CfeDfc ]^) = "R^)'^ 0) "^y 

Q.E.D. 

Corollary 3.5. For z ^ we have: 

(z m )'=mz m -\ m = -l,-2,.... 
Proof. Putting n = —771 and using Propositions 12.31 and I3T41 we obtain 



(z m Y =(— ] =-— (*")' — = -— nz"- 1 — = -77.: - = ///:' 

Q.E.D. 

Corollary 3.6. For an arbitrary constant c, we have: 

1 V 1 , 

r, z^c. 



c — z J (c — z)- 



Corollary 3.7. If quaternionic functions f and tp are M-differentiable at a 
point z° and ip 7^ in a neighborhood of z° , then the functions f ■ — and 
— ■ f are also M-differentiable at z° and we have: 

(f--)' ( z °) = /V) ■ ~hs - /(*°) "Ten ■ V'(*°) ■ "ToT 

V <pJ n* ) n z ) </?0 u ) 
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and 

(~ ■ /)' (*°) = -^oT ' <^°) ' -To^f^ + -Ten ■ /'(*°) • 

\ip J ip(z°) ip{z°) (^(z ) 

Proposition 3.8. Let a function f(z) be defined on some neighborhood of 
a point z° S H and let a function F(w) be defined on some neighborhood of 
the point w° — f(z). Assume that f is M- differ entiable at z° and that F 

is M-differentiable at w . If F'(w°) = ^2A k B k , then the composite Ff is 

k 
M-differentiable at z° and we have: 

{Ff)'{z ) = Y,Akf{z°)B k . 

k 

Proof. Let z be in the neighborhood of z . Put w — f(z). Then 
F(w) - F(w°) = J2 Ak( * w - w ^ Bk + ^( w °> w )> 

k 

f(z) f(z°) = J2 C & z )D j+U2 (z°,z), 

3 

and using these presentations, we calculate 

F(f{z)) F(f(z )) = J2 Mf(z) - f(z°))B k + wi(/( A /(*)) 

k 
k V j ' 

^^^■(z-^^ + oM+^it/tA/W)- 

k j 

But since 

M/( A /(*))! _ M/( A /(*))! \z-z°\ ;0) ^^ 
|z — z°\ \w — w°\ \w — w°\ 

we have 

(Ff)' (z°) = E E MCjDiBu = E A * ( E c ^) B * 

k j k ^ j ' 

= ^2A k f'(z°)B k . 



Q.E.D. 
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Specializing the proposition to the case where F(w) = w n and applying 
lUlJ) we get 



Corollary 3.9. // a function f is W-differentiable, then 

try - p- 1 ■ r + p- 2 •/'■/ + r- a • /' ■ f + ••• + /' • p- 1 - 

4 The H-derivative of the quaternion logarithm 
function 

A quaternion w is called the logarithm of a finite quaternion z ^ if z = e w , 
in which case we write w — In z . 

In order to define the H-derivative w' — (lnz)', we first note that the 
H-dcrivative of the left-hand side of the identity z = e lnz exits and is 1 by 
Proposition 12.31 Applying now Proposition 13.81 to the right-hand side and 
taking into account ()2.5[) . we get 



I , w w 2 \ , / 1 w w 2 \ , fl ,. 

l = (l + - + _ + ...j. w '+^+g f + ir + . ..].«/■« (4.1) 

1 W W 2 \ , 9 

3! + 4! + ^! + -"J- W - W +•••■ 
Thus, the H-derivative w' = (lnz)' satisfies Equality (|4.1I) . 

Remark 4.1. If vow 1 and w'w were equal, then we could write w ■ w', w 2 ■ 
w' , . . . instead of w' ■ w, w' ■ w 2 , . . . , and then Equality (|4.1j) would take the 
form 



/ ' w w 2 \ , ( w w 2 


■Yw' 


/ w 2 w 3 w A \ . 




= [l + w + — + — + ■■■ \-w' = e w -w' = 


■ e lnz •( 


So, we would obtain the classical formula 




(lnz)' = i, 

z 





(lnz)' = z- (lnz)'. 



that is well known in the case of a complex variable z. 
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